## Summary Statistics
## These results are reported in Table 1
rm(list=ls())
options(scipen=100,digits=10)
library(raster)
library(plyr)
library(dplyr)
library(lfe)
library(ggplot2)
library(stargazer)
load("data_final.RData")

#Change the variable names to get the summary statistics
mean(data3$SO2_MASS..tons.[which(data3$Year==2013)])
sd(data3$SO2_MASS..tons.[which(data3$Year==2013)])
mean(data3$SO2_MASS..tons.[which(data3$Year==2014)])
sd(data3$SO2_MASS..tons.[which(data3$Year==2014)])
mean(data3$SO2_MASS..tons.[which(data3$Year==2015)])
sd(data3$SO2_MASS..tons.[which(data3$Year==2015)])
mean(data3$SO2_MASS..tons.[which(data3$Year==2016)])
sd(data3$SO2_MASS..tons.[which(data3$Year==2016)])
mean(data3$SO2_MASS..tons.[which(data3$Year==2017)])
sd(data3$SO2_MASS..tons.[which(data3$Year==2017)])
mean(data3$SO2_MASS..tons.[which(data3$Year==2018)])
sd(data3$SO2_MASS..tons.[which(data3$Year==2018)])
mean(data3$SO2_MASS..tons.)
sd(data3$SO2_MASS..tons.)